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Motivated by the results of an experiment using atomic force microscopy performed by Gotsmann 
and Fuchs [Phys. Rev. Lett. 86, 2597 (2001)], where a strong energy loss due to the tip-sample 
interaction was measured, we investigate the potential implications of this energy loss channel to the 
quality factor of suspended micro- and nanoresonators. Because the observed tip-sample dissipation 
remains without a satisfactory theoretical explanation, two phenomenological models are proposed 
to generalize the experimental observations. A minimal phenomenological model simply extends 
for larger separations the range of validity of the power law found experimentally for the damping 
coefficient. A more elaborate phenomenological model assumes that the noncontact friction is a 
consequence of the Casimir force acting between the closely spaced surfaces. Both models provide 
quantitative results for the noncontact friction between any two objects which are then used to 
estimate the energy loss for suspended bar micro- and nanoresonators. Its is concluded that the 
energy loss due to the unknown mechanism has the potential to seriously restrict the quality factor 
of both micro- and nanoresonators. 



I. INTRODUCTION 

Noncontact friction, the energy loss due to the rela- 
tive motion of bodies moving without contact, has been 
receiving great attention in the context of atomic force 
microscopy (AFM)[l[. It has been investigated by its 
potential applications as a surface imaging technique, 
providing information about the surfaces not accessible 
by other techniques, as well as for the understanding of 
the basic physical mechanisms underlying the dissipation 
process. Most of the recent theoretical and experimen- 
tal efforts to understand this dissipation process were fo- 
cused on the noncontact friction at very short tip-sample 
separations, typically below 1 nm. While still a matter 
of debate, in this short distance (large force) regime, the 
evidences suggest that the fundamental physical mecha- 
nism leading to energy dissipation is force hysteresis 
Another important result regarding this subject was ob- 
tained by an earlier experiment which investigated the 
energy dissipation in AFM for comparatively larger tip- 
sample separations 0- More specifically, the dissipation 
between an aluminum coated silicon tip and a crystalline 
gold surface was measured and the results reported for 
tip-sample separation ranging from contact up to 4 nm, 
for the case of a tip with a radius i? = 35 ± 5 nm, and 0.5 
nm up to 7 nm for a tip with radius i? = 21 ± 2 nm. Un- 
der the assumption of viscous damping, the noncontact 
friction for separations larger than 1.5 nm was character- 
ized by a distance-dependent damping coefficient which 
scaled with tip-sample separation d as ^{d) cx d^^ . In 
spite of the attempts to explain the observed dissipa- 
tion, none of the known physical mechanisms leading to 
noncontact friction can explain the observed long range 
effects 0, 3- In fact, the measured dissipation is usually 
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orders of magnitude greater than any theoretical predic- 
tion . The existence of such noncontact friction for the 
aluminum-gold system rises questions about its potential 
implications for a variety of micro- and nanosystems or 
devices. The same large dissipation could exist for other 
material combinations, and its effects could extend over 
larger separations between material bodies. 

In the present work the dissipation due to an energy 
loss mechanism (ELM) described by two phenomenologi- 
cal models intended to generalize the results presented by 
Gotsmann and Fuchs[2] (GF) is investigated for rectan- 
gular cross section suspended beam micro- and nanores- 
onators. The loss of vibratory mechanical energy is one 
of the most important problems to be addressed into 
the current development of micro- and nanoresonators, 
specially for microelectromechanical (MEM) and nano- 
electromechanical (NEM) resonators. Several anelastic 
processes and clamping, for example, are known sources 
of energy dissipation 0, Q and have been investigated 
both theoretically and experimentally. However, it was 
only recently that noncontact friction was considered as a 
source of energy dissipation in such devices [1, Q- In gen- 
eral, the noncontact friction for suspended micro- and 
nanoresonators arises from the interaction of the res- 
onator with surrounding structures over vacuum or air 
gaps. Its relevance is evidenced, for instance, by the 
results presented in Ref. where, as a consequence of 
acoustic electromechanical energy loss, low quality fac- 
tors are predicted for realistic microresonators. 

Since no fundamental physical mechanism is known 
to explain the results of GF, the proposed ELM for the 
resonators investigated in this work is motivated by the 
analogy between the experimental setup of GF and a 
suspended resonator separated from nearby structures 
through nanogaps. In both cases we have surfaces in close 
proximity and in a periodic relative motion. Therefore, 
dissipation of the resonator mechanical energy can be ex- 
pected to be induced by the same mechanism present at 
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the tip-sample system in the typical configuration consid- 
ered here, where the resonator oscillates perpendicular 
to an underlaying plane surface (for example, the sub- 
strate or large electrodes). In principle, this analogy is 
strictly valid when the resonator and the nearby plane 
surface are made from the same materials used in the 
GF experimental tip-sample system. However, because 
the number o possible physical phenomena in any sim- 
ple system is limited, we can expect to observe the same 
physical mechanism inducing noncontact friction in anal- 
ogous systems involving other materials. In fact, at least 
for small tip-sample separations, energy dissipation of the 
same order of magnitude has been observed for other tip- 
sample materials combination involving metals, semicon- 
ductors and insulators T, "s*, 'o'l . It is this possibility that 
makes the present investigation more relevant, since the 
results and conclusions can be valid for other systems of 
practical interest. 

The article is organized as follows. In section |TT] we 
derive a suitable expression for the quality factor Q in 
terms of an areal damping coefficient. In section Hill two 
phenomenological models intended to explain and/or ex- 
tend the range of applicability of experimental results 
and their predictions for Q are presented. Further dis- 
cussions and our conclusions are presented in section HVl 



II. QUALITY FACTOR 



is characterized by a distance dependent areal damping 
coefficient r(d). In such a case, a resonator vibrating 
with a displacement Un{x,t) = u„(a;) sin(w„t) dissipates 
the power, 

P(t) = / T[d[x)\un[x,tf wdx 
Jo 

Ki r{do)wu!f^ Un{x)'^dx sin^iujnt), (1) 
Jo 

where the approximate expression is obtained by taking 
d{x) = do + Un{x) ~ do, because of the small amplitude 
of oscillation required for operation in the linear regime. 

Taking into account the orthonormality of the mode 
shapes Un{x) the average dissipated power results to be 

U^^T{do)wlLjlul, (2) 
leading to the quality factor 

The resulting quality factor depends on the geometry of 
the resonator only through the ratio h/l. It does not 
depend on the actual dimensions of the resonator, being 
the same for both micro- and nanoresonators. 



In this section we derive the expression for the qual- 
ity factor of a rectangular cross section suspended res- 
onator which loses energy due to a viscous damping 
force i^"««cotis _ ry(^(piy^ whcrc 7(c?) is the damping co- 
efficient. The resonator is considered to be homoge- 
neous, having length I, thickness h, and width w, each 
dimension being parallel to the x, y and z axis, respec- 
tively. It can be clamped or free at each end, result- 
ing in a cantilever, bridge or free resonator that vi- 
brates transversally in the z direction. It vibrates close 
to a plane surface contained in the xy plane and situ- 
ated at the average distance c?o- The mode shapes are 
those given by the Eulcr-Bernoulli beam theory, namely, 
Un{x) = uo{cosh{Knx/l) — cos{Knx/l) + [sinh(«;„a;//) — 
sin{Knx/l)]}, where n is the mode index (n = 1,2,3...), 
and Kn and Xn are determined accordingly for each 
boundary condition (Hj. 

The quality factor is defined as Q = 27rJ7„/AC/„ — 
i^nC/n/n, where Un denotes the total mechanical energy 
stored into the resonator, and A?7„ corresponds to the 
energy lost per cycle, which can be calculated from the 
time averaged energy loss 11. The vibrational energy 
for cantilevers, bridges and free resonators can be very 
well approximated by Un = where p de- 

notes the density and a;„ the mode frequency given by 
ujI = KiEh^/{12pl^), with E the Young modulus. The 
only missing ingredient for the calculation of Q is the av- 
erage energy loss. As we discuss further in the next sec- 
tion, we are going to assume that the noncontact friction 



III. PHENOMENOLOGICAL MODELS AND 
RESULTS 

A. Minimal phenomenological model 

1. Model 

In order to be able to use the experimental results for 
dissipation measured for a sphere-plane geometry (used 
to model the tip-sample system), for the approximately 
plane-plane geometry of the vibrating resonator and the 
nearby substrate, a fundamental assumption is that a 
physically meaningful areal damping coefficient T(d) can 
be defined. By physically meaningful it is meant that 
r(d) is the fundamental physical quantity determining 
the overall damping coefficient. 7(rf) is then obtained 
by performing a suitable integration limited to the over- 
lapping area of the two surfaces. We note that for sev- 
eral known noncontact friction mechanisms a physically 
meaningful T{d) can, in fact, be defined. That is the case 
for the van der Waals (vacuum) frictionQ or stochastic 
frictionflO]. However, for the electrostatic friction only 
j{d) can be defined due to the electric field dependence on 
the overall surface geometry, and similarly for the phonon 
emission mechanism, where the surface stress is, in gen- 
eral, a complex function of the applied forces 0, d, 0]- 

From the experimental data obtained by GF for the 
tip-sample damping coefficient "/{d) it is straightforward 
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to determine T{d). First, we note that for any two in- 
teracting surfaces the general relation between 7 and F 



lid) 



TdA, 



(4) 



where the integral is performed over the overlapping area, 
and d is a conveniently defined distance between the 
surfaces. In general, it is not straightforward to deter- 
mine F(d) from this relation, however, in the case of 
the tip-sample system we can resort to the Deryaguin 
approximation 
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We note that all the mathematical 
requirements for the use of the Deryaguin approximation 
in the evaluation of Eq. are satisfied, namely, F(c?) 
is a fast decreasing function of surface separation (based 
on fact that 7(d) (x d~^), and the tip has a radius of cur- 
vature R satisfying d, where d denotes the distance 
of closest sphere-plane separation. Using the Deryaguin 
approximation the approximate relation 



7(d) = 27ri? / T{z)dz 



(5) 



can be obtained, where z is the vertical coordinate per- 
pendicular to the plane surface. Deriving Eq. ([S]) with 
respect to d allow us to get. 



m = - 



1 d-f{d) 
2^ dd 



(6) 



It has to be noted that it is due to the particular sphere- 
plane geometry in the experiment of GF that we can eas- 
ily obtain the more fundamental physical quantity T(d) 
from the experimental results on ^(d). 

From the fitting of the experimental data on the aver- 
age dissipated energy (Figure 3 of Ref. as a function of 
the distance of closest tip-sample separation GF provided 
an expression for j{d) valid in the range 1.5 ^ c? ^ 7.0 
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where the approximate range of validity results from the 
uncertainties on the definition of the tip-sample distance 
inherent to the AFM technique. We can use this result 
and Eq. ^ to calculate T{d) in the range where the 
fit applies. However, the currently developed micro- and 
nanoresonators are made with gaps larger than approxi- 
mately 20 nm (see Ref. [H for the smallest gap produced 
so far) , and a model to extend the results of GF to larger 
distances becomes necessary in order to estimate the re- 
sulting quality factor for d > 7.0 nm. 

Based on the fact that T{d) or 7(0?) predicted by some 
of the known noncontact friction mechanisms also obey 
simple power laws on the variable d, it is reasonable to 
assume that the power law for ^(d) evidenced by the data 
fitting reflects the existence of an underlying physical 
mechanism that would be operating over larger distances. 
For this reason, we consider a minimal phenomenological 



model (MPM) which predicts T{d) by the substitution of 
7((i) into Eq. ([5]). The resulting damping coefficient is 



T{d) 



37" 1 
27ri?d4 



d^' 



(8) 



and is assumed to be valid for d > 1.5 nm. r° can be 
calculated taking into account the errors on the determi- 
nation of 7° = (8.011 j) X 10"^^ Nsm^ and i? = 21 ± 2 
nm. For an estimate of F^ we assume symmetrical er- 
rors for 7" equal to ±5 x 10"'^^ Nsm^ which results into 
ro = (1.8 ± 1.2) X 10-27 Nsm. 

We have just determined T{d) from measurements per- 
formed for a system comprised of two interacting sur- 
faces, one being crystalline gold and the other a not fully 
characterized aluminum surface 2] . While the results for 
y{d) can be considered as strictly valid for aluminum- 
gold systems, the analysis of the experimental data for 
a silicon-mica system performed in section IIVI indicates 
that the same ELM can exist for other material combi- 
nations. Assuming this point of view, we consider that 
T(d) obtained for the aluminum- gold system is only a 
reference value. In spite of the fact that no systematic 
experimental results obtained so far can exclude the role 
of the ELM described by the MPM on any other two 
materials combination, we can still be conservative and 
consider that this ELM is restricted to metallic systems. 
In this case, our results for Q would only be valid for 
whole metallic resonatorsfl^ or those, for instance, hav- 
ing metalized electrodes [14j. However, in what follows, 
no specific assumptions are made with regard to the exact 
structure of the resonator and the underlying substrate, 
and for the purpose of the mechanical modeling the res- 
onator is assumed to have a homogeneous and purely 
elastic structure. 



2. Results 

In Fig. [T] we present a range of Q values as a function 
of the gap d, for a doubly clamped polysilicon {E = 170 
GPa and p — 2.3 x 10"^ kg) resonator having a typical 
aspect ratio h/l = 0.1. We let r°, and consequently Q, 
vary over one order of magnitude around the reference 
value r" — 1.8 x 10^^^ Nsm. We have chosen polysil- 
icon as a reference material because it is the most fre- 
quently used structural material in micro- and nanode- 
vices. Values of Q in the same range are predicted for 
other structural materials because of the relatively small 
range of variation (within a factor of two) of the factor 
appearing in the Eq. (jSj for materials of relevance 
in micro- and nanofabrication including metals, semicon- 
ductors, and insulators. 

The predicted Q is rather small, specially below 20 
nm, even for the smallest value assumed for F'^. Such 
quality factors, well below 10^, are smaller than some of 
the lower values obtained experimentally for suspended 
bar nanoresonators not under the influence of nearby 
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FIG. 1. Log-linear plot of the quality factor Q predicted by 
the MPM as a function of the gap d for a doubly clamped 
polysilicon resonator with aspect ratio h/l = 0.1. The con- 
tinuous line corresponds to — 1.8 x 10~^^ Nsm. The up- 
per(lower) dashed lines are obtained by dividing(multiplying) 
F° by VW. 

surfaces [13I Therefore, this ELM has the poten- 

tial to significantly degrade Q if the currently fabricated 
nanoresonators are, for instance, set to operate as elec- 
tromechanical devices actuated electrostatically through 
nanometric gaps. It has also to be noted that, while the 
aspect ratio h/l = 0.1 is typical, much smaller aspect 
ratios can be easily found for practical devices, implying 
that the impact of this ELM on the degradation of Q 
can be even more significant. This ELM can also be rele- 
vant for MEM resonators, currently fabricated with gaps 
as small as about 30 nm[16], because the predicted Q is 
smaller than the normally found intrinsic quality factor 
of the order of 10"*. Another important result, which does 
not depend on the validity of the proposed MPM but, in- 
stead, depends solely on the validity of the assumption 
that there is a r{d) given by Eq. is the prediction of 
extremely small Q in the range supported by the exper- 
imental data d ^ 7 nm. This result indicates that non- 
contact friction, as measured by GF, can pose stringent 
limits on the quality factor of small gap nanoresonators 
or, equivalently, may set limits to the smallest possible 
practical gaps. 

B. Extended phenomenological model 

1. Model 

The results obtained based on the MPM demonstrate 
the potential relevance of the noncontact friction mech- 
anism unveiled by GF to both suspended micro- and 
nanoresonators. Therefore, a better understanding of 
this noncontact friction mechanism is desirable. In the 
search for a phenomenological model that could account 
for the experimental data presented by GF in a simple an 
coherent form we first note that the tip-sample force i^*" 
is essential in determining the friction resulting from dif- 



ferent mechanisms such as phonon emission[3|, stochastic 
friction^Oi], viscoelastic dissipation [l7|. and the simple 
model based on two-level systems of Ref. [H. In spite of 
the fact that none of those mechanism can explain the 
dissipation observed by GF they evidence the possible 
relation between i^*" and noncontact friction. The func- 
tional dependence of dissipation on the force can be quite 
involved as for the viscoelastic dissipation [l3|. It can also 
be straightforward as is the case for stochastic friction 
and the two-level systems mechanism proposed in Ref. 
[Tsl . In the first case it is predicted that 7(d) cx (VF**)^, 
while the last mechanism predicts 'y{d) oc F*''. 

Considering that in the experiment of GF the role 
of electrostatic force was minimized living the Casimir 
force as the only relevant force at large tip-sample 
separations, we further analyzed the data presented by 
GF and established a possible relation between the tip- 
sample force and the damping coefficient. We note that 
j{d) given in Eq. 1^ depends on d as d~^ in a distance 
range where the Casimir force reduces to the van der 
Waals force of the form F""^^ = -HR/(&d^), where H 
denotes the Hamaker constant [ll|, and whose gradient 
ypts ^ dF*'/dd = HR/{M^) also goes as d'^. This 
fact suggests the possibility that 7(d) be proportional to 
the force gradient. Based on this possibility we formu- 
late the following hypothesis for the relation between the 
tip-sample force and the damping coefficient 



where C is a phenomenological constant. 

Unfortunately, the data provided by GF are presented 
for experiments involving AFM tips with different radius 
and do not form a coherent data set that allow us to test 
our hypothesis over a broader range of probed tip-sample 
distances. The data for the energy dissipation presented 
m Fig. 3 of Ref. i in the range 0.5-7.0 nm , and fitted in 
the range 1.5-7.0 nm leading to Eq. Q, is not accompa- 
nied by the corresponding data on the force that could be 
used to calculate the energy dissipation predicted based 
on the hypothesis Eq. However, another data set 

for the measured tip-sample force and the corresponding 
^{dY^'P (obtained from a special numerical fitting of the 
experimental data on the dissipated power |2|), referring 
to measurements performed using a tip with i? = 35 ± 5 
nm, was presented by GF in Fig. 1. This figure pro- 
vides the curves that allow for a test of the hypothesis 
we formulated in the range of short distances below 1.5 
nm where the force increases exponentially. The data 
obtained from the curves using a high resolution auto- 
mated digitizer is not sufficiently precise for an analysis 
for d ^ 1.5 nm, where the force starts its transition to 
pure van der Waals force. From the extracted data points 
for the experimental force we numerically calculate the 
force gradient, the data being smoothed by the method 
of moving averages. By comparison with the curve for 
^{dy^^ we determined the value of the phenomenolog- 
ical constant as C^'-^" = (l.l ± 0.7) x 10"^ s, where 
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d (nm) 

FIG. 2. Log-linear plot of the experimental 7(d) as extracted 
from Fig. 1 of GF (continuous line), and the theoretical 7(d) 
calculated from the experimental data for the tip-sample force 
and Eq. ((9]) (dashed line). 



the error was estimated from the combined error in the 
calculation of C = 7((i)'^^P/VF*'', with the uncertainties 
on -f{dyP and the force being taken as 40%[2,]. The 
experimental and theoretical curves can be seen in Fig. 
[21 where the results are presented in the range 0.5-1.5 
nm. For d < 0.5 nm 7(0?)'^^^ keeps growing exponentially 
while the theoretical 7(rf) tends to decrease. The good 
superposition of the experimental and theoretical curves 
serves to corroborate the hypothesis represented by Eq. 
([9]), which seems to be valid for d > 0.5 nm. Below 0.5 
nm some other noncontact friction mechanism that does 
not complies with our hypothesis can dominate the dis- 
sipation. 

Motivated by the seeming validity of the hypothesis 
represented by Eq. ([9]) , we propose a second phenomeno- 
logical model that not only explains the measured dissi- 
pation in the range 0.5 < d < 7.0 nm but is also as- 
sumed to be valid for larger distances. The extended 
phenomenological model (EPM) is then based upon the 
previous assumption that ^{d) derives from the more fun- 
damental physical quantity T{d) with the addition of the 
new hypothesis. Now, taking into account that Eq. ([9]), 
so far, is considered strictly valid for a tip-sample system, 
we can apply the Deryaguin approximation to rewrite 
both 7(d) and F'** in such a way that 



7(d) = C 



Qpts 

~dd 



r°° d r°° 

2ttR / T{z)dz = 2ttRC— / 
Jd od Jd 



'■{z)dz, (10) 



where F^'^{z) denotes the pressure between two flat sur- 
faces separated by a distance z. Consequently, T{d) and 
F'^" (d) must be related by 



r(d) = C 



dF^^id) 
dd ' 



(11) 



Before we proceed, we note that this fundamental rela- 
tion, obtained by performing an approximate calculation 



to the particular tip-sample system, could be derived 
more rigorously if the hypothesis expressed in Eq. (jH]) 
was generalized by stating that the damping coefficient 
is directly proportional to the force gradient for any two 
interacting surfaces. This is a reasonable but not neces- 
sary generalization, and we do not discuss it further. 



2. Results 

To calculate Q predicted by the EPM we need to know 
C, and determine F** (d) for our specific system. Because 
we are interested in systems with gaps larger than a few 
nanometers, F'^'^d) is dominated by the Casimir force as 
predicted by the Lifshitz theorv [l9l[20j . This force results 
from the quantum fluctuations of the vacuum electromag- 
netic field, and is a function of the optical properties of 
the surfaces. It can be conveniently expressed as 



F-(d)=,7(d) 



hcTT 

2A0d* 



vid)F^id), 



(12) 



where 77(d) corresponds to the correction factor to the 
Casimir force between two perfectly conducting surfaces 
F'-^(d). In the present analysis, this correction factor 
takes into account only the effects of the finite conduc- 
tivity of the Au and Al surfaces, and was calculated us- 
ing the Lifshitz theory as described, for instance, in Ref. 
I20I . The optical data for both metals is that provided 
by Palik[2l|. At small separations of a few nanome- 
ters ry(d) oc d, and the Casimir force reduces to the van 
der Waals force. Having obtained Hamaker constant, 
jjAi-Au 3 g ^ io~^^ J, we can calculate the force be- 
tween two plane surfaces and determine C in an alter- 
native manner from the same data used to evaluate Q 
predicted by the MPM. Substituting F*"(d) = F""^^ {d) 
in Eq. and taking the experimental result for 7(d) 
from Eq. (O, valid for the short distance limit, we ob- 



tain C 



Al-Au 



(3.0 



-1-2.6 



1.8)X10" 



s, where the errors derive 



from the uncertainties on 7" and R. While the central 
value differs significantly from the previously obtained 



value of C 



Al — Au 



(1.1 ± 0.7) X 10"® s, both results can 



be reconciled if we take into account the estimates for 
the errors. 

Because the value C^'--^" = (3.0+^J) x 10"® s was de- 
termined from the same 7(d) used to calculate r(d) pre- 
dicted by the MPM, an adequate comparison between the 
predictions of both MPM and EPM requires the adop- 
tion of this value for C and not the value obtained from 
other data set. For the sake of comparison Q predicted 
by the EPM was calculated for the same resonator con- 
sidered for the MPM using r(d) given in Eq. ((TT|) with 
F'"'{d) calculated numerically from Lifshitz theory, and 
using as a reference value (7^'^-^" — 3.0 x 10~® s. The 
results gained from the EPM are compared with those 
predicted by the MPM in Fig. [3] As expected from the 
faster decrease of the force predicted by the Lifshitz the- 
ory when compared to the van der Waals force, the dissi- 
pation predicted by the EPM becomes smaller than that 
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FIG. 3. Log-linear plot of the quality factor Q predicted by 
the EPM as a function of the gap d for a doubly clamped 
polysilicon resonator with aspect ratio h/l = 0.1. The con- 
tinuous line corresponds to C"^'--*" = 3.0 x 10~* s. The up- 
per(lower) dashed lines are obtained by dividing(multiplying) 
qAi-Au ^jg rpj^g dot-dashed line presents Q predicted 
by the MPM and corresponds to the continuous line in Fig. 

m 



for the MPM for larger distances resulting into higher 
quality factors. It can be seen that because the Casimir 
force reduces to the van der Waals force varying as at 
short separations the predictions for Q coincide for small 
d. 

At this point, it is worth to compare our results with 
those obtained from another phenomenological model of 
noncontact dissipation proposed earlier by Dedkov(22|. 
The model is proposed in an attempt to get an unified 
explanation for the experimental results for both paral- 
lel and transversal relative motion between tip and sam- 
ple. This model assumes that the damping coefficient is 
directly proportional to the tip-sample force instead of 
the force gradient, that means ^{d) — V^^F*''{d), where 
y is a phenomenological velocity to be determined from 
the experimental data. The model was tested using the 
data on tip-sample force and damping coefficient for the 
Al-Au system investigated by GF and for the silicon- 
mica system investigated in Ref. i9|, in the case of relative 
transversal motion. The analysis was restricted to short 
tip-sample separations, smaller than 2 nm, and "f{d) pre- 
dicted based on the experimental data on the tip-sample 
force agreed very well with the experimental data on j{d) 
with an adequate choice of the phenomenological veloc- 
ity. As we discuss further in the next section, the predic- 
tions of the EPM for 'y{d) using the experimental data 
on the force for the silicon-mica system are also in good 
agreement with the measured ^{d). The reason both phe- 
nomenological models can account for the experimental 
data can be attributed to the fact that the measured 
and 7((i) are exponentially varying functions of 
d in the range 0.5 < d < 2.0 nm. The EPM never- 
theless can also account for the experimental results for 
d > 2 nm, where the phenomenological model proposed 
by Dedkov predicts incorrectly that 7(d) cx d~^ for the 




0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

d (nm) 



FIG. 4. Data points and best fit curves (/3 — 2.35 nm~^) 
for IF*'* I (•/continuous line) and 7 (o/dashed line) for the 
sihcon-mica system. 

Al-Au system. 



IV. DISCUSSION AND CONCLUSIONS 

The results presented so far for Q based on the mini- 
mal and extended phenomenological models highlight the 
potential relevance of an yet unknown energy loss mech- 
anism on the performance of micro and nanoresonators. 
However, the determination of the phenomenological con- 
stant and the corroboration of the EPM relied on exper- 
imental results for an Al-Au system, therefore restrict- 
ing the generalization for other systems of practical in- 
terest, as mentioned at the Introduction. Fortunately, 
the experimental results on force and dissipation for a 
semiconductor-insulator system is also available that fur- 
ther corroborate the EPM. In Ref. l9|the tip-sample force 
and dissipation were measured for a silicon tip (i? ~ 10 
nm) vibrating close to a mica surface, the data presented 
in the Fig. 7 of Ref. ^ for F*" (d) and 7(d) being accurate 
for a comparison between theory and experiment for tip- 
sample separations limited to less than approximately 2 
nm. The data points for the force were extracted using 
an automated digitizer and are presented in Fig. |4l The 
experimental data for the damping coefficient was con- 
veniently extracted from the work of Dedkov [2 2], where 
it is reproduced, and is also presented in Fig. |4l As 
for the Al-Au system in the short range limit, the best 
fit for both force and damping coefficient for d in the 
range 0.5-2.0 nm were obtained with exponentially vary- 
ing functions of the form Q;exp(— /3d), where a and f5 are 
the fitting parameters. While the best fit value of (3 for 
F**(d) and 7(d) differ by about 10%, adequate goodness 
of fit and predictions differing from the best fit values by 
less than the estimated experimental errors are obtained 
by the assumption, consistent with the EPM, of an in- 
termediate value of /3 = 2.35 nm~^ for both F*''(d) and 
7(d). The resulting best fit curves are plotted in Fig. SI 

The phenomenological constant obtained from the 
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substitution of the best fit functions in Eq. ([9]) was 

(jSi-niica ^ Q 94 ^ ^ j.gg^i|. comparable to that 

for the Al-Au system for the short-range data. More 
precisely, performing the same analysis performed above 
for the silicon-mica system the best fit to the short- 
range data for the Al-Au system results in an inverse 
decay length /3 = 2.2 nm~^, similar to that for the 
silicon-mica system. Furthermore, the phenomenolog- 
ical constant obtained by this procedure results to be 
IjAi-Au ^ 3 ^ iQ-?, ^ j.gg^i|. compatible with that 

derived directly from the experimental data by the curve 
matching in section UlI B II 

From this discussion we can conclude that the EPM is 
corroborated by experiments involving very distinct ma- 
terials combination. Additionally, the phenomenological 
constant that characterizes the yet unknown underlying 
physical mechanism has similar values for both materi- 
als combination, indicating that the fundamental process 
leading to energy dissipation could be the same irrespec- 
tive to the materials involved. However, further experi- 
mental investigations of noncontact friction are required, 
specially for larger separations, because several questions 
are still open regarding, for instance, the role of spe- 
cific physical conditions, chemical composition, rough- 
ness, temperature and geometry of the surfaces. 

The experimental results should provide a guidance for 
more fruitful theoretical investigations that goes beyond 
the phenomenological models. In fact, we note that the 
results from a set of five experiments already point to- 
ward a fundamental role of surface imperfections [1, 2, ^, 
IHj nil on energy dissipation. All five experiments in- 
volved the interaction of an imperfect tip surface with, 
nominally perfect, atomically flat surfaces. In the three 
experiments where the tip vibrates perpendicular to the 
sample surface (T], d, Q and, consequently, is subject to a 
time varying force damping is several orders of mag nitude 
larger than that for the two experiments [23l, 'm| where 
the tip vibrates parallel to the sample and is subject to a 
constant force. Because in both configurations the sam- 
ple surface is always subject to a time varying force it 
is reasonable to conclude that the main contribution to 
the very distinct damping levels comes from energy dis- 
sipation due to the time varying force acting at the tip 
surface. In fact, the experiment by Kisiel et a/.[2j], has 
already elucidated that in the case of parallel oscillation 
the main contribution to noncontact friction comes from 
electronic friction, a mechanism that can not explain the 
dissipation for vertical oscillations. The role of the tip 
surface is further corroborated by the fact that dissipa- 



tion has shown no dependence on sample temperature in 
the experiment reported in Ref. [ll, a result that alone has 
led the authors of this reference to consider that the ma- 
jor contribution (to dissipation) stems from mechanisms 
within the tip[l]. It has to be noted that in this last work 
experimental results for the tip-sample force and dissi- 
pation are also presented for silicon tips on crystalline 
NaCl(OOl) surface. However, the range of useful data for 
analysis ranges from 0.5 nm up to only approximatelly 
1.2 nm and leads to results that are inconclusive with 
regard to the relation between F**'(d) and 7(d). Such 
a result is expected if the tip radius becomes too small 
compared to i? ~ 20 — 35 nm in the experiment by GF, 
because for small tip radius the role of the atoms at the 
tip apex can dominate the overall tip-sample force and 
dissipation[25] at the very small tip-sample separations 
where they are experimentally accessible. In this case the 
prevailing mechanism of energy dissipation seems to be 
force hysteresis as evidenced by atomistic simulations f25j 
and experimental results [26l[27|. 

Finally, from the results for the quality factors pre- 
dicted based on both the MPM and the EPM, we con- 
clude that the yet unknown ELM first observed by GF 
can pose stringent limits on two of the most relevant op- 
erational parameters of both micro- and nanoresonators, 
the quality factor and the minimum gap. For instance, 
because low mechanical impedances and high Q are usu- 
ally sought in MEM and NEM resonators for RF appli- 
cations, a trade off between these two operational pa- 
rameters may be a consequence of this ELM, for the 
impedance and Q increase rapidly with the gap. Fur- 
ther experimental investigations aimed at measuring dis- 
sipation at larger separations up to tens of nanometers 
should provide the pieces of information required to test 
the two phenomenological models, and for gaining a bet- 
ter understanding of the long range noncontact friction 
and its consequences for micro- and nanodevices. 
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